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We solve the RSOS(p) models on the light-cone lattice with fixed boundary 
conditions by disentangling the type II representations of SU(2) q , at q = e l7V / p , 
from the full SOS spectrum obtained through Algebraic Bethe Ansatz. The rule 
which realizes the quantum group reduction to the RSOS states is that there must 
not be singular roots in the solutions of the Bethe Ansatz equations describing 
the states with quantum spin J < (p — l)/2. By studying how this rule is active 
on the particle states, we are able to give a microscopic derivation of the lattice 
S— matrix of the massive kinks. The correspondence between the light-cone Six- 
Vertex model and the Sine-Gordon field theory implies that the continuum limit 
of the RSOS(jo + 1) model is to be identified with the p— restricted Sine-Gordon 
field theory. 
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1. Introduction 



It is known that the Six- Vertex (6V) model, in the so-called light-cone for- 
mulation and with periodic boundary conditions (p.b.c), yields the Sine-Gordon 
massive field theory in an appropriate scaling limit [1]. Hence the light-cone 6V 
model can be regarded as an exactly integrable lattice (minkowskian) regularization 
of the SG model. 

Recently, the hidden invariance of the SG model under the quantum group 
SU(2) q was exhibited [2]. In the first part of this work [3], henceforth denoted as 
paper I, we showed how this SU {2) q symmetry can be enforced also in the regular- 
ized lattice version, by choosing suitable fixed boundary conditions (f.b.c.) on the 
light-cone 6V model. In this formulation, the quantum group generators J z , J+ 
and J_ appearing in the limit of infinite spectral parameter of the Yang-Baxter al- 
gebra generators, can be explicitly shown [3] to commute with the transfer matrix. 
Moreover, the algebraic Bethe Ansatz (BA) remains applicable [4] and one can 
verify that the f.b.c. BA states are highest weight states of SU(2) q [3,5]. It follows 
that the Solid-On-Solid (SOS) states are, by definition, just the f.b.c. BA states. 
The 6V Hilbert space includes the whole SU(2) q multiplets and follow by repeat- 
edly applying the lowering operator J_ to the highest weight BA states. In paper I 
we also generalized the light-cone approach to this f.b.c. contest, constructed the 
massive scaling limit of the 6V and SOS models and derived the higher-level BA 
equations describing the internal space degeneracy of the physical excitations. 

In this second part we complete our work as follows. After a brief summary 
of the results of paper I (section 2), we present in section 3 and 4 the derivation 
of the factorized S— matrices on the lattice, i.e. still in the presence of the UV 
cutoff. This derivation is based on the "renormalization" of the BA Equations, 
which consists in removing the infinitely many roots describing the ground state. 
What is left is once again a f.b.c. BA structure involving the lattice rapidities of 
the physical excitations (the particles of the model) and the roots of the higher- 
level BAE obtained in paper I. The explicit form of two-body S— matrix for the 
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6V model and the SOS model can be extracted in a precise way (eft. eq. (3.11)) 
from this higher-level BA structure. In the massive scaling limit these lattice 
scattering amplitudes become the relativistic 5 1 — matrices of the SG model (or 
Massive Thirring model) and of the continuum SOS model. Let us remark that 
the SOS S— matrix, although closely related to the 6V and SG S— matrices from 
the analytical point of view, is conceptually different. It describes the scattering 
of kinks interpolating between renormalized local vacua labelled by integers. This 
kink S— matrix is most conveniently expressed in the so-called face language (see 
eq. (4.3)). 

In section 5 we investigate the f.b.c. BAE (eq. (2.4)) when the quantum group 
deformation parameter q is a root of unity, say q p = ±1, with p some integer 
larger than 2 (the case p = 2 being trivial). In this case it is known that RSOS(p) 
model can be introduced by restricting to the finite set (1, 2, . . . , p — 1) the allowed 
values of the local height variables of the SOS model [6]. This is equivalent to the 
projection of the full SOS Hilbert space (which is formed by the highest weight 
states of SU(2) q ) to the subspace spanned by the type II representations [7], that 
is those representations which remain irreducible when q becomes a root of unity. 
Our results on this matters, in the BA context, can be summarized as follows: 

a) Only when q is a root of unity, the f.b.c. BAE admit singular roots (that is 
vanishing z— roots or diverging v— roots, in the notations of sec.l). 

b) When q tend to a root of unity, say q p = ±1, the BAE solutions can be 
divided into regular and singular solutions, having, respectively, no singular 
roots or some singular roots. Regular solutions correspond to irreducible type 
II representations. Singular solutions with r singular roots correspond to 
the reducible and generally indecomposable type I representations obtained 
by mixing two standard SU(2) q irreps of spin J and J + r (we recall that 
J = N/2 — M, where N is the spatial size of the lattice and M the number 
of BA roots). Then necessarily r < p. 
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c) The r singular roots zi, Z2, ■ ■ ■ , z r vanish with fixed ratios 

Zj=J~ 1 z 1 , uj = e 27Ti/r , l<j<r (1.1) 

In terms of the more traditional hyperbolic parametrization, with v— roots 
related by Vj = — \ log^j to the z— roots, the singular roots form an asymp- 
totic string-like configuration. They have a common diverging real part and 
are separated by n/r in the imaginary direction. 

So we see that the RSOS eigenstates are easily singled out from the full set 
of BA eigenstates of the 6V or SOS transfer matrix. One must retain all and 
only those BAE solutions with M > (N — p + l)/2 (which correspond to states 
with J < (p — l)/2) and with all M z— roots different from zero. This provides 
therefore an exact, explicit and quite simple solution for the RSOS model on the 
lattice (with suitable boundary conditions, as we shall later see). In particular the 
ground state of the 6V, SOS and RSOS models is the same f.b.c. BA state (in the 
thermodynamic limit N — > oo at fixed lattice spacing). It is the unique SU(2) q 
singlet with all real positive roots and no holes. The local height configuration 
which, loosely speaking, dominates this ground state can be depicted as a sequence 
of bare kinks jumping back and forth between neighboring bare vacua (see fig. 1). 
In the massive scaling limit proper of the light-cone approach [1], this ground state 
becomes the physical vacuum of the SG model as well as of and all the Restricted 
SG field theories. 

In the discussion ending sec. 5, we argue that the kink S— matrix for the 
excitations of the RSOS models follows indeed by restriction from that of the SOS 
model. In the scaling limit it is to be identified with the relativistic S— matrix 
of the Restricted SG models. All these field-theoretical 5 1 — matrices are naturally 
related to the Boltzmann weights of the respective lattice models. This is because 
the higher-level BAE are identical in form to the "bare" BAE, apart from the 
renormalization of the anisotropy parameter 7 (related to q by q = e* 7 ) 



and the replacement of the rapidity cutoff ±0 with the suitably scaled rapidities 
9j of the physical excitations 

±e ^ Jh_ (1.3) 

7T — 7 

In particular, since 7 = 71 /p for the RSOS(p) model, eq. (f .2) yields p — > p — 1 . 
This shows that the renormalized local vacua £ run from 1 to p — 2 when the 
bare local heights i n run from 1 to p — 1. The higher-level BA structure of the 
light-cone 6V model (or lattice regularized SG model) thus provides a microscopic 
derivation of the bootstrap construction of ref. [2] and explains why the S— matrix 
of the p— restricted SG field theory has the same functional form of the Boltzmann 
weights of the lattice RSOS(j» + 1) model. Moreover, the well-known correspon- 
dence between the critical RSOS(p) models and the Minimal CFT Models M p 
imply the natural identification of the massive p— restricted SG model with a com- 
pletely massive relevant perturbation of M p . This is generally recognized as the 
perturbation induced by the primary operator 0^3 with negative coupling. 

Two detailed examples of the BA realization of the quantun group reduction 
to the RSOS models are presented in section 6. We considered the simplest cases 
p = 3 and p = 4. The RSOS (3) model is a trivial statistical system with only 
one state, since all the local heigths are fixed once we choose, for example, the 
boundary condition £q = 1. In our quantum group covariant f.b.c. construction 
this corresponds to the existence of one and only one type II state when 7 = ir/3. 
We then obtain the following purely mathematical result: for any N > 2 and real 
w = exp(-2Q) the set of BAE 



Zjiv-e™/ 3 zj -we™/ 3 \ N _ y? ] Zj - z k e 2%i ^ Zj z k -e 2%i / 3 i < <N 

ZjWe 7ri/3 _ 1 Zje 7ri/3 _ W J 11 Zje 27vi/3 _ Zk ZjZke 2iri/3 _ 1 ' ~ J ~ 

(1.4) 

admit one and only one solution with non-zero roots within the unit disk \z\ < 1. 
In addition, these roots are all real and positive. 
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The RS0S(4) model can be exactly mapped into an anisotropic Ising model as 
showed in eqs. (6.1-5). In our case the horizontal and vertical Ising couplings turn 
out to be 0— dependent complex numbers. For even N, the Ising spins are fixed on 
both space boundaries. For odd N the spins are fixed on the left and free to vary 
on the right. We analyzed the BAE for the RSOS(4) model in some detail. In the 
thermodynamic limit the ground state, as already stated, is common to the 6V and 
SOS moodels. The elementary excitations correspond to the presence of holes in 
the sea of real roots charactering the ground state. Each holes describes a physical 
particle or kink and may be accompanied by complex roots. In sec. 5 we argue 
that in the RSOS(4) case a state with v holes necessarily contains \v/2\ two-strings 
those position is entirely fixed by the holes. Notice that here the number v of holes 
can be odd even for even N. This is not the case for the 6V or SOS models, where 
v is always even for N even. What happens is that when 7 — > (7r/4) _ the largest 
real v— root diverges in the J = 1 states of the 6V and SOS models. Therefore, 
these states get mixed with J = 2 states into type I representations and do not 
belong to the RSOS(4) Hilbert space. The RSOS(4) states with J = 1 are obtained 
by choosing the largest quantum integer I N/2-1 = N/2 (see the Appendix). There 
is no root associated to N/2 + 1. It follows that these states, from the 6V and SOS 
viewpoint, have a cutoff dependent term in the energy equal to ir/a, where a is the 
lattice spacing (loosely speaking, one could say that "there is a hole at infinity"). 
They are removed from the physical SG spectrum in the continuum limit. We are 
thus led to propose as RSOS(4) hamiltonian, for even N 

H RSOS(4) = Hsosil = tt/4) - a^nJ (1.5) 

where H$os is given by eqs. (2.9), (2.1), (2.16-20) and J = or 1. In this way, 
the particle content of the light-cone RSOS(4) and the corresponding S- matrix 
turn out to coincide with the results of the bootstrap-like approach of ref. [2]. 
Eq. (1.5) defines, in the scaling limit, the hamiltonian of the (p = 3)— restricted 
SG model. Notice, in this respect, that the higher-level BAE (2.12) and (3.10) 
completely determine the physical states in terms of renormalized parameters. 
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To summarize, the picture we get from the BA solution of the f.b.c. 6V, SOS 
and RSOS(p) lattice models is as follows. Performing the scaling limit whithin the 
light-cone approach, these lattice models yields respectively: the SG model (or 
Massive Thirring Model), a truncated SG and the (p — l)restricted SG models. 
For the SG model we have essentially nothing to add to the existing literature, 
apart from the explicit unveiling, at the regularized microscopic level, of its SU (2) q 
invariance and for a better derivation of the S— matrix. The truncated SG follows 
by keeping only the highest weight states with respect to SU (2) q , that is the kernel 
of the raising operator J + . Finally we showed that the RSOS(p) lattice models with 
trigonometric weights yield in the scaling limit proper of the light-cone approach 
the (p — l)restricted SG field theories formulated at the bootstrap level in ref. [2]. 



2. Summary of the Quantum Group invariant 
formulation of the light-cone Six-vertex model 

In paper I we formulated the 6V model on a diagonal lattice with special 
fixed boundary conditions. This conditions are such that the dynamics of the 
model enjoys explicitly quantum group invariance under SU(2) q . The unit-time 
evolution operator U on the lattice with N diagonal links, can be written 



U = -R12-R34 • • • RN-l-e,N~e9N9l 1 #23-^45 • • • -RiV-2+e,iV-l+e 



(2.1) 



where e = [1 — (— l) N ]/2, gj is the matrix exp0<r z acting nontrivially only on the 
two-dimensional vector space attached to the j th link, and Rjk is the 6V R-matrix 
R(2Q) acting nontrivially only on the tensor product of the j th and k th vector 
spaces. In the conventions adopted in this paper, this R-matrix reads, in terms of 
a generic spectral parameter 0, 



R(0) 



(\ 0\ 

c b 

b c 

\fl 1/ 



c =c(0,7) 



sinh# 



sinh(i7 — 0) 

sinh i'-f 
sinh(i7 — 0) 



(2.2) 
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where the anisotropy parameter 7 (we may assume < 7 < 7r) is related to the 
quantum group deformation q by q = exp(i7). 

The R-matrix (2.2) is unitary for real 9. Hence the evolution operator U is 
unitary too, for real 6, up to the boundary effects due to the term g^g^ 1 in eq. 
(2.1). The exact diagonalization of U, to be described below, shows that all its 
eigenvalues are unimodular, so that there exists a similarity transformation map- 
ping U to an explicitly unitary operator. On the other hand, it is natural to expect 
that in the thermodynamic limit (N — > 00) different boundary conditions become 
equivalent (the model has a finite mass gap). This suggest that the above men- 
tioned similarity transformation reduces to the identity as iV — > 00. We conclude 
therefore that the light-cone 6V model with fixed b.c. described by the evolution 
operator U of eq. (2.1) is another good, integrability-preserving regularization of 
the SG model. Its advantage over the more conventional setup with periodic b.c. 
is that U is explicitly SU(2) q invariant even for finite N (see paper I for details). 

The Algebraic Bethe Ansatz leading to the exact diagonalization of U (together 
with the whole family t(9) of row-to-row transfer matrices with inhomogeneities 
±0) is described in detail in paper I. The main results are as follows: 

1. The eigenvectors of U are written 

v 2 , ...,v M ) = B^Bfa) ■ ■ ■ B(v M )n (2.3) 

where Q is the ferromagnetic reference state and the operators B(vj) create 
exact multi-spin-wave states with J z = N/2 — M labeled by the unordered 
set of distinct "bare rapidities" v\, . . . ,vm, < M < N/2. 

2. The exponentials Zj = exp(—2vj) must satisfy the set of algebraic Bethe 
Ansatz Equations 

/ \ N M 2 2 

f ZjW-q Zj-wq \ yv Zj_ - z k q z Zj z k - q z 

\zjwq - 1 zjq-wj £1 Zjq 2 _ Zk ZjZkq 2 -\ 
where w = exp(— 26). As explained in paper I, the search for z— roots of 
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these BAE can be restricted in the interior of the unit circle. 
3. The eigenvalue of U on the BA state (2.3) reads 



, ¥t sinh(9 - Vj + zt/2) sinh(6 + Vj + i 7 /2) 

u(vi,V2,...,v M ) = 11 . , =2 ■ /0 v ■ r / n , : TTyT (2.5) 

f-^ sinh (6 — — -ry/2) smh(9 + -Uj — ry/2) 



4. The BA states are highest weight states of SU(2) q , in the sense that they 
are annihilated by the raising operator J + . Hence to each BA state, that is 
to each solution of the BAE (2.4), there corresponds an entire multiplet of 
states: for irrational values of 7/77, i.e. for q not a root of unity, this multiplet 
has dimension 2J + 1, with J = J z , and all other lower weight states with 
J z < J are obtained from the J z = J BA state by repeated application of 
the lowering operator J_. The highest weight J is related to the eigenvalue 
of the q— Casimir C q on each multiplet: 

C q = Sin(J+1 ^ 7 SinJ7 (2.6) 
sin 7 

Let us recall that C q appears in the limit of infinite spectral parameter of the 
transfer matrix t(9), 

f(+oo) = 2 cos 7 -4C g sin 2 7 (2.7) 

while, in the same limit 

B(+oo) = (l-q 2 )J- . (2.8) 

According to the general light-cone construction, the physical lattice hamilto- 
nian H is defined in terms of U in the natural way 

U = e- iaH (2.9) 

where a is the lattice spacing (the same in both space and time directions). Then 
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the BA states have energy 



M 

E = a^^eoivj) 



(2.10) 



where 



eo(^) = + 2 arctan 



( 



cosh 20 cos 7 — cosh 2v 



sinh 20 sin 7 



) 



(2.11) 



Notice that eo(v) is smooth and negative definite for real v. The specific choice 
of the logarithmic branch in passing from eq. (2.5) to eqs. (2.10), (2.11) is dictated 
by the requirement that eo(f) should correspond to the negative energy branch of 
the spectrum of a single spin wave over the ferromagnetic state Q. 

The physical ground state and the particle-like excitations are obtained by 
filling the interacting sea of negative energy states. This sea is described by a set 
of real v— roots of the BAE with no "holes" (in the Appendix we give a detailed 
exposition on the treatment of the BAE, also to clarify some rather subtle matters). 
Excitations correspond to solutions of the BAE which necessarily contain holes and 
possibly complex v— roots. The crucial point is that in the limit iV — > 00 only the 
number of real roots of the sea grows like N, while the number of holes and complex 
roots stays finite to guarantee a finite energy above the ground state energy. Hence 
these solutions of the BAE can be described by densities p(v) of real roots plus a 
finite number of parameter associated to the positions x\,X2, ■ ■ -x v of the v holes 
and to the location of the complex v— roots. In paper I we showed that, in the 
limit N 00, these complex roots can be organized into "2-strings", "wide pairs" 
and "quartets" , described by a set of parameters Xi, X25 • • • Xm with M = u/2 — J. 
Most importantly, these parameters depend on the positions of the holes through 
the higher-level, or renormalized, BAE 



h=i 



n 



jjWh - q Q - w h q 
ijW h q - 1 ijq - w h 




- ikf - q 2 

s/r - i jil-q 2 - 1 



(2.12) 
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where 

= e" 2Q » , w h = e- 2aXh , q = e* = e ia ~< , a = (2.13) 

These new BAE involve only a finite number of parameters and are formally iden- 
tical to the "bare" BAE (2.4), with the holes acting as sources in the place of the 
alternating rapidities ±0. Moreover, from (2.12) one reads out the renormalization 
of the quantum group deformation parameter 

q^q i.e. 7^7 = (2.14) 

7T — 7 

As we shall see later on, this renormalization has a nice physical interpretation for 
the SOS and RSOS models related to the 6V model by the vertex-face correspon- 
dence. 

For completeness, we close this summary with a brief description of vertex-face 
correspondence just mentioned. On the faces of the diagonal lattice one introduces 
positive integer-valued local height variables £ n (t), t E 2L, according to (see fig. 2) 

£ (t) = l; h(t)=2 

(2.15) 

e n+1 (t)=£ n (t)±l, n = l,2,...N-l 

The configurations {£i(t), . . . ^iv(O) at an Y fixed discrete time t are in one-to- 
one correspondence with the SU(2) q multiplets of the 6V model of dimension 
£ n (t) = 2J+ 1 (notice that £jv(t) can be chosen to be time-independent thanks to 
SU(2) q invariance). Now the matrix elements of U between these highest weight 
states define the unit-time evolution in the face language. However, in order to 
write such matrix elements in a a simple factorized form, one needs to pass first to 
a different vertex representation by means of the similarity transformation 

U = GUG- 1 (2.16) 

with 

n 1/2 -1/2 £-1/2 , . 

G = 9i 9 2 ■■■9n (2.17) 
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Then one explicitly finds 



?(„ 4 ■■■ 41 ^ Ko, 4, . . . <?jv} = w^i'w^ . . . w' N _ 1+e w 2 W4 . . . w N _ e (2.18) 



where the W n are the face weights 



w n = w(e n - h e n+1 \e n ,e' n ) ■ w' n = ( 2 -i9) 



and 



W^(i,j|r, s) = 5 rs + 5 ij b{2Q, 1 ) 



sm7r sin 7s 
sin 72 sin 7,7 



1/2 



(2.20) 



3. Higher-level Bethe Ansatz and S-matrix 

The holes in the sea of real BAE roots are the particles of the light-cone 6V 
model. They are SU(2) q doublets and can be present in even (odd) number for iV 
even (odd). Consider first the even iV sector, setting iV = 2N' . The energy of a 
BA state with an even number v of holes located at x\, ■ ■ ■ x v can be calculated 
to be, in the N' — > oo limit, 

Ef cosh it^c / '"y "\ 
e(x h ) + 0(a- 1 N~ 1 ) , e(x)=2arctan — - — (3.1) 
\ smh 7T0 / 7 / 

where Eq is the energy of the ground state, that is the state with no holes. Eq is 
of order and is explicitly given in the appendix along with some detail on the 
derivation of eq. (3.1). Now suppose odd, with = 2A^' — 1. To compare this 
situation with the previous one, we need to slightly dilate the lattice spacing a to 
a' = 2N'a/(2N' — 1), in order to keep constant the physical size L = Na of the 
system. Then eq. (3.1) remains perfectly valid, as N' — > oo, also for the case of 
odd, with the same ground state energy Eq. The only difference is that now v 
is odd. Hence, altogether, we obtain that the number of holes can be arbitrary 
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(unlike in the treatment with periodic b.c.) and that the total energy, relative to the 
ground state and in the L — > oo limit, is the uncorrelated sum of the energy of each 
single hole, independently of the complex pair structure of the corresponding BAE 
solution. The BA state with one hole has J = N/2 - M = (2N 1 - 1)/2 - (TV' - 1) = 
1/2 and is therefore a SU(2) q doublet. The 2 V polarizations of a state with v 
holes are obtained by considering all solutions of the higher-level BAE (2.12) with 
< M < v/2. We see in this way that the holes can be consistently interpreted 
as particles. 

As a lattice system the light-cone 6V model is not critical. The (dimensionless) 
mass gap is the minimum value of the positive definite e(x), the energy of a single 
hole, that is 

e(0) = 2 arctan ( , 1 ) (3.2) 



sinh 7r6/7 

This gap vanishes in the limit — > oo. Hence the continuum limit a — > can be 
reached provided at the same time 6 — > oo in such a way that the physical mass 

m = a- 1 e(0) ~ Aa^e^^ (3.3) 

stays constant [1]. In the same limit we obtain the relativistic expression 

a~ 1 e(x) — > mcosli7ra/7 (3.4) 

so that nx/'-f is naturally interpreted as the physical rapidity 9 of the hole. This 
is consistent also in our fixed b.c. framework, provided the limit L — > oo is taken 
before the continuum limit. Indeed the total momentum becomes a conserved 
quantity in the infinite volume limit and its eigenvalues can be expressed in terms 
of the BA v— roots exactly as in the periodic b.c. formulation. Then one finds 
that the momentum of a single particle takes the required form msinh# in the 
continuum limit. 
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The above analysis shows that a relativistic particle spectrum appears in the 
a — > limit above the antiferromagnetic ground state. For 7 > n/2 one can 
show also that bound states appear, associated to appropriate strings of complex 
roots, just as in the periodic b.c setup. In the sequel we shall anyway restrict our 
analysis to the repulsive 7 < tt/2 region, where the only particles are the holes, to 
be identified with the solitons of the SG model. Of course at this moment, since 
the infinite volume limit is already implicit, the particles are in their asymptotic, 
free states: the rapidities 9h = KXh/l ma Y assume arbitrary continuus values and 
the total excitation energy is the sum of each particle energy and does not depend 
on the internal state of the particles. The situation changes if we consider L very 
large but finite, since in this case the hole parameters still quantized 

through the "bare" BAE (2.4). Indeed, by definition the holes are real distinct 
numbers satisfying 

ZN(xh',v\,V2,...,v M ) = ; h = l,...,v (3.5) 

where Z^(x ; v\, . . . , vm) is the "counting function" defined in the Appendix and 
the positive integers I\, . . . , l v are all distinct from the integers I\, . . . , I r labelling 
the r (r < M) real v— roots of the BAE. For iV very large, J = N/2 — M finite 
and x < (7/vr) IniV, the counting function can be approximated as 

Z N (x;vi, . . . ,vm) = Z OQ (x) + N~ 1 F(x]Xi,. . . ,xm]Xi,- ■ ■ ,Xm) + ° { n ~ 2 ) ( 3 -6) 

where Z^ix) is the ground state counting function at the thermodynamic limit 

r - r , _ / sinli7r:r/7\ . . 

Zoo (x) = 2 arctan — ' 3.7 
\ v cosh7rB/7/ 

and 

V 

F(x ;xi,..., x M ;Xi, Xm) = ~i log JJ S (^(x - x h ))S (%(x + x h )) 

h=\ 

M ( 3 - 8 ) 

t-i- sinh a(x — Xj + *7/2) sinho;(x + Xj + *7/2) 
— i log I I 

sinho;(a; — Xj — «7/2) sinho;(x + Xj — «7/2) 
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In the last expression, the numbers Xj are the roots of the higher level BAE (2.12), 
(2.13), while Sq(9) coincides with the soliton-soliton scattering amplitude of the 
SG model 



Sq(9) = exp % 



dk sinh(7r/27 — l)k sin k9/7r 
k sinh(7r/c/27) cosh/c/2 



(3.9) 



under the standard identification 7/71 = 1 — (3 2 /8tt. 

Combining eqs. (3.5) and (3.6), and taking the continuum limit one obtains 
the "higher level" expression 



M 



exp (-imLsmh9 h ) = J] S (9 h - 6 n )S (6 h + 9 n ) J{ 



n=l 



jjWh - q £jq ~ Wh 



(3.10) 



where, according to (2.13), £j = e~ 2aXj and = e~ 2 ^ 6h /' K . Together with the 
higher level BAE (2.12), this last equation provides the exact Bethe ansatz di- 
agonalization of the commuting family formed by the v renormalized one-soliton 
evolution operators (see fig. 3) 



Uh =Sh,h-l(0h - ) • • • Sh,l{6h - 81) 9h(—0h) Sh,l(9fi + &l) ■ ■ ■ Sh^(9 h + 9„)x 
9h{0h) S h; „(9 h - 0„) . . . S h;h+ i(9 h - 6h+i) 

(3.11) 

where gh{9) = exp 9a? and S(9) is the complete 4x4 SG soliton S— matrix in the 
repulsive regime 7 < 7r/2 (for brevity we set 9 = j9/(tt — 7)): 



S(9) = S (9) 



(\ 0\ 

6 c 

c 6 

\fl 1/ 



6=6(0,7) 



c =c(6,i) 



sinh 9 



sinh (27 — 9) 

sinh 
sinh (27 — 9) 



(3.12) 



The operators Uh are the values at 9 = 9h of the fully inhomogeneous Sklyanin- 
type transfer matrix T(9 ; 9\, . . . , 9 V ) constructed with S(9 — 9h) as local vertices. 
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Let us stress that the higher-level Bethe Ansatz structure just described follows 
directly, after specification of the BAE solution corresponding to the ground state 
and without any other assumptions, from the "bare" BA structure of the light- 
cone 6V model. In particular, this provides a derivation "from first principles" of 
the SG S— matrix. We want to remark that S(9) is also the exact S— matrix for 
the elementary excitations of the 6V model on the infinite lattice: it is a bona- 
fide lattice S— matrix. Of course, in this case 9 = 9\ — 62 is the difference of 
lattice rapidities, which are related to the energy and momentum of the scattering 
particles through the lattice uniformization 



Thus, we see that the continuum limit only changes these energy-momentum re- 
lations to the standard relativistic form, (Ej, Pj) = m(cosh9j, sinhfy), without 
affecting the S— matrix as a function of 9. 

4. SOS and RSOS reductions and kink interpretation 

Under the vertex-face correspondence previously described, the light-cone 6V 
model is mapped into the SOS model. At the level of the Hilbert space, this 
corresponds to the restriction to the highest weight states of SU(2) q : each (2 J + 
1)— dimensional multiplet of spin J (for generic, non-rational values of 7/71) is 
regarded as a single state of the SOS model. In particular, in such a state the local 
height variables £ n have well defined boundary values £q — 1 and £n = 2J + 1. The 
ground state of the 6V model, which is a SU(2) q singlet, is also the ground state 




(3.13) 



implying the lattice dispersion relation 



cosaEj/2 = tanh(7r©/7) cosaPj/2 



(3.14) 
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of the SOS model: it has £q — £n — 1 and is "dominated", in the thermodinamic 
limit N — > oo by the see-saw configuration depicted in fig. 1. Our boundary 
conditions allow for only one such "ground state dominating" configuration, with 
£n = (3 — (— ) n )/2, while periodic b.c. on the SOS variables £ n would allow any 
possibility: £ n = (21 + 1 ± (— ) n )/2, with £ any positive integer. 

Consider now a BA state with one hole. Then N is necessarily odd while 
J = 1/2 and £n = 2 J+l = 2. From the SOS point of view this state is "dominated" 
by height configurations of the type depicted in fig. 4. Upon "renormalization" , we 
can replace the ground state and one hole configurations with the smoothed ones 
of fig. 5. The hole corresponds, in configuration space, to a kink of the SOS model 
which interpolates between two neighboring vacuum states: the vacuum on the 
left of the kink has £ n = (3 — (— ) n )/2, corresponding to a constant "renormalized" 
height £ n = 1, while the vacuum on the right has £ n = (5 + (— ) n )/2, corresponding 
to £ n = 2. Of course, many kink configurations like that depicted in fig. 4 are to 
be combined into a standing wave (a plane wave in the infinite volume limit which 
turns the one hole BA state into an eigenstate of momentum). Let us observe 
that, by expressions like "dominating configuration", we do not mean that, e.g. 
the ground state, becomes an eigenstate of £ n as N — > oo. We expect local height 
fluctuations to be present even in the thermodynamic limit or, in other words, that 
the ground state remains a superpositions of different height configurations. The 
identification of a dominating configuration is made possible by the integrability 
of the model, which guarantees the existence of the higher-level BA. In turns, the 
higher-level BA allows us to consistently interpret the ground state or the one hole 
state as in figs. 1, 4 and 5, since the presence of fluctuations only renormalizes in a 
trivial way the scattering of physical excitations relative to the bare, or microscopic, 
R— matrix (2.2), as evident from eq. (3.12). Therefore, we can reinterpret at the 
renormalized level the vertex-face correspondence: the holes, that is the solitons of 
the SG model, are SU(2) q doublets acting as SOS kinks that increase or decrease 
by 1 the renormalized heights £ n . The higher-level BA makes sure that the total 
number of internal states of v kinks interpolating between £ = 1 and £ = 2 J is just 
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the number of highest weight states of spin J in the tensor product of v doublets 

^ = ©"(/- 1) ; k ^^- J <"> 

This parallels exactly the original BA, which provides the 

dN(J) = (m) - (m - l) ; *=W2]-J < 4 - 2 > 

highest wight states of N doublets with total spin J. 

Comparing eqs. (2.2), (2.20) and (3.12), we can directly write down the 
S— matrix of the SOS kinks: 

5(C;j=5o(fl)|*„ + 5y 6(^,7) 

It defines the two-body scattering as follows. The first ingoing kink interpolates 
between the local vacuum with £ — i and that with £ = r (r — i ± 1), while the 
second interpolates between £ = r and £ — j (j — r ± 1). The outgoing kinks are 
interpolating between £ — i and £ = s (s — i ± l,i.e, s = r, r ± 2) and between 
£ = s and I = j (j = s ± 1). 

When 7 = n/p for p = 3,4,..., the SOS models can be restricted to the 
RSOS(p) models, by imposing £ n < p. In the vertex language of the 6V model 
this corresponds to the Hilbert space reduction to the subspace formed by the 
so-called type II representations of SU(2) q with q p = — 1 [7]. In the next section 
we shall describe in detail how the restriction takes place in our BA framework. 
Here we simply observe that the local height restriction, when combined with the 
finite renormalization 7 — > 7 of eq. (2.14), provides a strong support for the kink 
interpretation presented above. Indeed, if 7 = n/p, then 7 = 7r/(p — 1) and the 
renormalized heights £ n can take the values 1, 2, ... ,p — 2. This appears now ob- 
vious, since each constant configuration of £ n corresponds to an £ n — configuration 



siii^r shits 
sin ji sin7j 



(4.3) 
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oscillating two neighboring values. Moreover under the standard identification of 
the critical RSOS(p) models with the minimal CFT series M p , we see that each 
light-cone RSOS(p) model has kink excitations whose S— matrix (4.3) is propor- 
tional to the (complex) microscopic Boltzmann wights of the RSOS(p— 1) model, 
under the replacement of26 = — (— 6)= rapidity difference of light-cone right 
and left movers, with (p— 1) times 6 = 61—62 =rapidity difference of physical parti- 
cles. This is exactly the pattern found by bootstrap techniques [2] for the minimal 
model M p perturbed by the primary operator 0i.3 (with negative coupling). 



5. BA roots when q is a root of unity 
and Quantum Group reduction. 

As previously explained to each M roots solution of the BAE there corresponds 
a highest weight state of the quantum group SU{2) q with spin J = N/2 — M. It 
is well known that when q is a root of unity, say q p = ±1, then {J+) p = {J-) p = 
and the representations of SU(2) q divide into two very different types. Type I 
representations are reducible and generally indecomposable. They can be described 
as pairwise mixings of standard irreps (that is the irreducible representations for q 
not a root of unity) with spin J and J' such that 

\J-J'\<p, J+J'=p-1 (modp) (5.1) 

Notice that the sum of q— dimensions for this pair of reps vanishes. 

Type II representations are all the others. They are still fully irreducible and 
structured just like the usual SU{2) irreps. Since (J±) p = 0, type II representations 
have necessarily dimension smaller than p, that is 

J < ^ (5-2) 

We shall now show how the BAE reflect these peculiar properties of the SU(2) q 
representations for q a root of unity. First of all let us stress that when q is not a 
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root of unity [i.e.^jn is irrational), then the BAE cannot possess v— roots at (real) 
infinity. Indeed, since B(oo) is proportional to J_ (eq. (2.8)), a root at infinity 
means that the corresponding BA state is obtained by the action of the lowering 
operator J_ on some other state with higher spin projection J z . But for q not 
a root of unity, the BA states have J = J z and therefore cannot be obtained by 
applying J_ on any other state. On the other hand, if we assume that one z— root, 
say zi, of the BAE (2.4) lay at the origin, i.e.Revi = +00, then eqs. (2.4) for 
j = 1 imply 

g 4 ( J+1 ) = 1 (5.3) 

That is, q must be a root of unity. It is now crucial to observe that the remaining 
equations for the non-zero roots (those labelled by j = 2, 3, . . . , M) are precisely 
the BAE for M — l unknowns. This invariance property holds only for the quantum 
group covariant, fixed boundary conditions BAE. It does not hold for the p.b.c. 
BAE where v— roots at infinity twist the remaining equations for finite roots. This 
twisting reflect the fact that the corresponding periodic row-to-row tranfer matrix 
is not quantum group invariant but gets twisted under SU(2) q transformations [8]. 
Thus, when q is a root of unity, BA states (2.3) with one v — root at infinity take 
the form 

= 00, v 2 , . . . , v M ) = (1 - q 2 ) J-*(«2, V3, . . . , vm) (5-4) 

Let us recall that the BA state on the l.h.s. is annihilated by J + for any q (including 
q a root of unity) [3]. Hence eq. (5.4) represents the mixing of two reps with spin 
J and J' = J + 1 into a type I representation. Indeed, applying the mixing rule 
(5.1) into the necessary condition (5.3) for one root at infinity, yields an identity 
as required: 

1 = g4(J+l) = q 2(J+J'+l) = {qP) 2n = ^ (55) 

where n is a suitable (J— dependent) positive integer. 

We can generalize this analysis to any number of vanishing BA z— roots. Let 
us identify the roots going to the origin with z\, Z2, ■ ■ ■ , z r , 1 < r < M. The BAE 
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for the remaining non-zero roots z r +i, . . . , zm take the standard form (2.4) valid 
for a BA state formed by M — r B— operators, which has therefore spin J + r. The 
BAE for the vanishing roots take the form 

g 4J+2(r+l) = p TT ZjQ 2 ~ *k X < J < r (5 6) 

k& 

where the definition of Fj is understood in the limit of vanishing z\, . . . , z r . We can 
obtain complete agreement with the quantum group mixing rules (5.1) by setting 
r < p and making the natural choice 

Fj = l, 1 < j < r (5.7) 

Indeed the presence of r vanishing z— roots imply the mixing of two representations 
with spin J and J' = J + r, so that 

g 4J+2(r+l) = g 2(J+J'+l) = {g) 2n = j (5g) 

with n an integer depending on J and J'. Eq. (5.7) ha a very simple solution for 
the limiting behaviour of the vanishing roots. We find that if 

zj = ujJ" 1 2i , 1 < j < r (5.9) 

with a> r = 1, then identically 

- ain[- 7 + 7 r(Jb-j)/r] = 
J sm[ 7 + 7r(A;-j)/r] 

for any value of 7. That eq. (5.10) should hold for generic values of g = e* 7 is 
necessary since we are studying the approach of the roots to the origin when 7/7? 
tends to a rational number. We would like to remark that the limiting behaviour of 
the ratios of the vanishing BA z— roots depend only on their number r and neither 
on the specific rational value of •y/n nor on the other non-vanishing roots. 
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In the v— plane, the roots vi, V2, ■ ■ • , v r go to infinity as a r— string with spacing 
m/r. For example, when r = 2 we have a pair of complex roots with imaginary 
parts tending to ±27r/4 as their real parts simultaneously diverge. For r = 3 there 
will be a limiting 3-string with a real member and a complex pair at ±iir/3. For 
these two examples, we have performed an explicit numerical test when 7 — > (n/3)~ 
and 7 — > (7r/4) ~, respectively, finding perfect agreement with the picture proposed 
here. We shall present more detail on the special cases in the next section. Let us 
anticipate here that the numerical results suggest the following conjecture on the 
behaviour of the singular roots: \zj\ ps 0(e 1//r ) and lm[vj — Vj(e = 0)] ps 0(e 1 ^ r ) 
where 6 = 7 — n/(r + 1). 

When q p = ±1 and the SU(2) q reps divide into type I and type II, it is possible 
to perform the consistent RSOS(p) reduction of the Hilbert space. This consists 
in keeping, among all states annihilated by J+ which form the SOS subspace, only 
the type II states. It is known that this corresponds to SOS configurations with 
local heights £ n restricted to the set 1,2, ...,p— 1 [7]. In our BA framework, 
therefore, the RSOS reduction is obtained by retaining only those states with spin 
J < (p-l)/2, that is with M > (N -p+l)/2, such that all BA roots zi,z 2 , ...,z r 
are non-zero. This constitues the general and simple prescription to select the 
RSOS(p) subspace of eigenstates of the evolution operator when q p = ±1. 

We have considered up to here the RSOS reduction at the microscopic level, 
that is for the BA states (2.3) described by the bare BAE (2.4). It should be 
clear, however, that the analysis of the singular BAE solutions when q is a root 
of unity holds equally well for the higher-level BAE (2.12), due to their structural 
identity with the bare BAE. The crucial problem is whether a quantum group 
reduction carried out the higher level would be equivalent to that described above 
in the "bare" framework. This equivalence can actually be established as follows. 
Suppose that 7/7? is irrational, but as close as we like to a specific rational value. 
Then all solutions of the bare BAE are regular and can be correctly analyzed, in 
the limit N — > 00, using the density description for the real roots. As described 
in paper I, this yields the higher-level BAE (2.12), those roots are in a precise 
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correspondance with the complex roots of the bare BAE. When 7 tends to ir/p, 
then 7 tends to ir/ (p — 1) and singular solutions of the bare BAE with two or more 
singular roots are in one-to-one correspondance with singular solutions of the 
higher-level BAE. Thus in this case quantum group reduction and renormalization 
of the BAE commute. The only potentially troublesome cases are those of singular 
solutions with only one real singular root. Indeed the density description of the 
v— roots cannot account, by construction, for real roots at infinity, and no sign of 
the type I nature of the corresponding BA state would show up at the higher level. 
Notice that spin J BA states with a single singular root are mixing with spin J + 1 
states, so that, by the quantum group mixing rule (5.1), necessarily J = p/2 — 1. 
On the other hand, a direct application of the constraint (5.2) at the higher level, 
that is with the replacement p — > p — p — 1, would overdo the job, by incorrectly 
ruling out all BA states with J = p/2 — 1 = (p — l)/2. To be definite, consider N 
even. Then the type II BA states with J = p/2 — 1, which ar e superpositions of 
SOS configurations with l n < p — I = 2J + I = l^, have an effective, higher-level 
spin J = (£n — l)/2 = J— 1/2. This higher-level spin does satisfy (5.2). The type 
I states are those in which £ n somewhere exceeds p — 1. In particular, the simplest 
change on a type II configuration, turning it into type I, is to flip the oscillating 
£ n close to the right wall (as shown in fig. 6). £n and hence the spin J are left 
unchanged, but clearly In increases by one, causing J to violate the bound (5.2). 
In other words, RSOS(p)— acceptable BA states with J = p/2 — 1 have one kink 
less than the corresponding SOS states. From the detailed analysis of the p = 4 
case, to be discussed below, it appears that the removal of one kink corresponds to 
giving infinite rapidity to the hole representing that kink in the higher-level BAE. 

Thank to the possibility of performing the restriction directly at the renormal- 
ized level, the kink S— matrix of the RSOS(p) model follows by direct restriction 
on the SOS S— matrix given by eq. (4.3). Namely, one must set 7 = ir/p, that is 
7 = 7r/(p — 1), and consider all indices as running from 1 to p — 2. 
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6. The models RSOS(3) and RSOS(4) 



In order to clarify matters about the BA RSOS reduction discussed in the pre- 
vious section, we present here the more details about the two simplest examples, 
when 7 = 7r/3 and 7 = 7r/4. For these two cases the RSOS reductions correspond, 
respectively, to a trivial one-state model and to the Ising model (at zero external 
field and non-critical temperature). Let us recall that the light-cone approach 
yields in the continuum limit massive field theories with the same internal sym- 
metry of the corresponding critical regimes. Therefore the RSOS (4) reduction of 
the light-cone 6V model coincides with the 7L2— preserving perturbation of the 
c = 1/2 minimal model (which is obtained upon quantum group restriction of the 
critical 6V model). 

6.1. The case p = 3 

From the RSOS viewpoint, the case 7 = n/3 is particulary simple. Since p = 3 
the local hight variables £ n can assume only the values 1 and 2. Then the SOS 
adjacency rule \£ — 1'\ = 1 implies that the global configuration is completely de- 
termined once the height of any given site is chosen. In our light-cone formulation 
the first height on the left is frozen to the value 1 (see eq. (2.15) and fig. 2), so 
that the restricted Hilbert space contains only one state: an SU(2) q singlet when 
N is even and the spin up component of a J = 1/2 doublet when N is odd. 

For 7 = 7r/3 the BAE still have many solutions which reproduce the full SOS 
Hilbert space. Indeed, from the SOS point of view, nothing particular happens 
when 7 — * 7r/3. However, at this precise value of the anisptropy, only one BAE 
solution corresponds to the unique type II representation: for N even (odd) it 
contains N/2 (N/2 — 1/2) non-zero roots in the z— plane. All other BAE solutions 
with M = [N/2] contain at least one vanishing z— root and correspond to type I 
SOS states. Moreover, there unique RSOS state is the ground state of the SOS 
model and is formed by real positive roots labelled by consecutive quantum integers 
(no holes). We reach therefore the following rather non-trivial conclusion: the 
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complicated system of algebraic equations (2.4) admit, for q = ±1 and w real, 
one and only one solution with M = [N/2] non-zero roots within the unit circle. 
In addition, these roots are real and positive. For few specific choices of N we also 
verified this picture numerically. 

6.2. The case p = 4 

The local height variables can assume now the three values £=1,2, 3. However, 
each configuration can be decomposed into two sub-configurations laying on the 
two sublattices formed by even and odd faces, respectively. On one of the two, the 
SOS adjacency rule \£ — £'\ = 1 freezes the local heights to take the constant value 
£ = 2. Then on the other sub-lattice we are left with two possibilities, £ = 1 and 3. 
Moreover, the interaction round-a-site of the original model reduces in this way 
to a nearest-neighbor interaction in the vertical and horizontal directions. The 
framework is that of the Ising model. 

To obtain the standard Ising formulation, we can set 

a n {t)=£ 2n {t)-2 (6.1) 

where the numbering of the lattice faces can be read from fig. 2. The fixed b.c. on 
the £ n now correspond to 

<7 (t) = -l, a R (t) = (-l) J - 1 (N = 2R) 

(6.2) 

<7 (*) = -l, a R (t) = ±l (N = 2R+1) 

where J = 0, 1 for N even, due to the bound (5.2). For N odd we must consider 
only the possibility J = 1/2, which implies £jy = 2, since the line of half-plaquettes 
on the extreme right belong to the frozen sublattice. Then £n~i is left free to 
fluctuate between 1 and 3, leading to free b.c. on a. The matrix elements of the 
unit time evolution operator U can be calculated from the explicit form of the SOS 
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weights (2.19). They can be written in the "lagrangian" form 



t+i 



(a' ,a' 1 ,...a' R \U\ao,a 1 ,... a R \ = e iL ^ (6.3) 



where 

R-l R-l 

L{t) = P v Vn(t)a n (t + l)-f3 h J2 <7n(t)<T n +i(t) + const (6.4) 

n=l n=0 

and 

V = J(tt + 2?lntanh29) , (3 h = arctantanh20 (6.5) 

Alternatively, standard simple manipulations allow to rewrite U explicitly in terms 
of Pauli matrices 

jj =e -il3 h H 2 e ~if3 h Hi 

n—1 n—Q 

In either cases, one sees that the complex Boltzmann weights formally belong to 
the critical line 

sin 2/^ sin 2/^ = 1 (6.7) 

Of course, this follows from the original definition of the light-cone 6V model in 
terms of complex trigonometric Boltzmann weights, which, under the replacement 
— > iO, would correspond to the critical standard 6V model. Nevertheless, just as 
for the vertex model, also in this "light-cone" Ising model, a massive field theory 
can be constructed in the Re 6 — > +oo limit. 

In the BA diagonalization of U, we must restrict ourselves to the BAE solutions 
with M = R or M = R — 1 for even iV and M = R for odd N. For sufficently large 
even N, the ground state has R real roots with consecutive quantum integers. 
Actually, as already stated above, this is a general fact valid for all RSOS(p) 
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models. Namely, the infinite volume ground state of th SG model, of the SOS 
model and of all its restrictions RSOS(p) is the same f.b.c. BA state. It is the 
unique SU(2) q singlet with all real positive roots and no holes. It is described in 
the thermodynamic limit (N — > oo, a fixed) by the density of roots given in eq. 
(5.7) of paper I. 

Excited states with J = have an even number of holes and a certain number 
of complex roots such that no z— root lays at the origin. For instance a two- 
particle state contains two holes (holes are naturally identified with the particles) 
and a two-string with imaginary parts ±[7r/8+ corrections exponentially small in 
N] laying between the two holes. This state (that is this precise choice of quantum 
integers) is just the two-particle state of the SG or SOS model, with 7 fixed to the 
precise value 7r/4. We have explicitly checked, by numerically solving the BAE for 
various values of N, that indeed all z— roots stay away from the origin as 7 crosses 
7r/4 while the quantum integers are kept fixed to the two-hole configuration. Now 
consider a state with four holes. Apart from the multiplicity of the rapidity phase 
space, there are two distinct type of such states: in v— space one contains two 
two-strings, while the other contains one sigle wide pair, that is a complex pair 
with imaginary part larger than 7. This situation holds for generic values of 7/7? 
and simply reflects the fact that the holes 

are SU(2) q doublets. The crucial point is that, as 7 reaches 7r/4 (from below), 
the wide pair moves towards infinity, while its real part gets closer and closer to 
the (diverging) value of the largest real root and the imaginary part approaches 
7r/3. This picture is confirmed by a careful numerical study of the BAE and is in 
perfect agreement with the general picture presented in sec. 5. In addition, there 
exists numerical evidence that the real parts diverge like log(7r/4 — 7)~ 1//6 while 
the imaginary part of the wide pair goes to tt/3 like (it/4 — 7) 1 / 3 . When 7 > 7r/4, 
the v — root corresponding to the largest quantum integer has the largest real part, 
but is no longer real, having an imaginary part equal to tt/2. As 7 — > (7r/4) + , 
this real part again diverges together with real part of the wide pair. Hence the 
four-hole state with a wide pair is type I when 7 = it/ 4. On the other hand, the 
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four-hole state with two two-string is type II, since all its v— roots stay finite. 

From the study of the two- and four-hole states, we are led to the following 
general conjecture: in the BA framework, the RSOS(4) J = states are all and 
only the states with v — 2k holes and k two-strings. Then, in the higher-level 
BAE (2.12), we recognize this state as that corresponding to the unique solution 
with k real roots. In other words, this BA state is completely determined once 
the location of the holes (that is the rapidities of the physical particles) is given. 

Consider now the states with J = 1. In this sector, the lowest energy type II 
state contains exactly one hole. This is an unusual situation for BA systems on 
lattices with even N, where holes are always treated in pairs. As long as 7 7^ 7r/4, 
the same is true in our f.b.c. BA: the lowest energy state with J = 1 contains two 
holes, in agreement with the interpretation of holes as SG solitons with quantum 
spin 1/2. For 7 < 7r/4 all roots are real. For 7 > 7r/4 the root fjv/2-1 corresponding 
to the largest quantum integer In/2-1 = ^/2 + 1 aquires an imaginary part 71/ 2 
(see the appendix for details). But when 7 = 7r/4 then ReVjv/2-1 — +00, and the 
J = 1 two-hole state is mixed with some J = 2 state into a type I representation. 
It does not belong to the RSOS(4) hilbert space. This picture can be easily verified 
numerically. 

To prevent the largest root v N/2-1 from diverging, it is sufficent to consider 
a J = 1 state with only one hole and In/2-1 — N/2. From the point of view of 
the SG or SOS models this one-hole state has a cutoff-dependent energy which 
diverges as a -1 in the continuum limit. If x\ is the position of the single hole, the 
energy relative to the ground state reads 

E - Eq = a _1 e(xi) + 7ra" 1 (6.8) 

where the renormalized energy function e(x) is given in eq. (3.1). Notice that 
this result coincides with the limit X2 — > 00 of a two-hole state (eft. eq. (3.1)). 
Strictly speaking, however, the density method leading to (3.1) has no justification 
when "one hole is at infinity". We trust eq. (6.8) nonetheless because it passes all 
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our numerical checks. In the continuum limit a — > 0, — > oo, 4a _1 e~ 7I " e> / 7 = m 
(fixed), this one-hole state is removed from the physical spectrum, as required 
from the SG point of view: holes are spin 1/2 solitons while here J = 1. 

The preceding discussion easily extends to generic states in the J = 1 sector 
containing an odd number of holes and a given set of complex pairs. The infinite- 
volume energy of these multiparticle states is given by 

V 

E = Eo + a' 1 ^ e(x h ) + vra" 1 (6.9) 
h=\ 

with the same divergent constant ixa^ 1 appearing irrespective of the physical con- 
tent of the state. Our conjecture for the J = states naturally extends to these 
J = 1 states: if there are u = 2k+l holes, then the v— roots are all finite, implying 
that the state is type II, provided there are also k two-strings. Then the higher- 
level BAE imply that these states are completely defined by the hole rapidities. To 
retain these type II J = 1 BA states in the RSOS(4) model, an extra J— dependent 
subtraction is necessary to get rid of the divergent constant. Namely, for N even 
and J = or 1, we set 

h rsos(4) = H S os(l = tt/4) - 7TGT 1 J (6.10) 

where H sos is given by eqs. (2.1), (2.9), (2.16-20). 

Alltogether, we see that the BA picture for the excitations of the RSOS(4) 
model is fully consistent with the kink interpretation for the holes. Indeed J = 
corresponds to even Dirichlet b.c. for the Ising field, while J = 1 corresponds to 
odd Dirichlet b.c. (eft. eqs. (6.2)). In the Ising model the kink description is valid 
below the critical temperature, with the disorder field as natural interpolator for 
the kinks. In this case the S— matrix must be —1 (it is +1 when the asymptotic 
particles are the free massive Majorana fermions) and this is exactly what follows 
from eq. (4.3) upon setting 7 = n/3 and restricting all indices to run from 1 to 
p-2 = 2. 
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The analisys of the BA spectrum for odd N does not contain real new features. 
Now J is fixed to 1/2 and the b.c. on the Ising field are of mixed fixed-free type, 
as shown in the second of eqs. (6.2). The lowest energy state corresponds to the 
BAE solution formed by real roots and one single hole as close as possible to the 
the v— origin. In other words, the quantum integers are given by Ij — j + 1, for 
j — 1, 2, . . . , (N — l)/2. By letting this hole to move away along the positive real 
axis, we reconstruct the energy spectrum of a one-particle state. Of course, to 
compare this state to the global ground state, which contains no holes, a judicious 
choice of the odd value of N is required. If in the ground state N = 2R, then 
the new state is indeed an excited state if we choose N = 2R — 1, rather than 
N = 2R + 1, because of the antiferromagnetic nature of the interaction. 



APPENDIX 

In this appendix, we present for completeness a (rather non-standard) treat- 
ment of the BAE (2.4). As explained in paper I, it is convenient to first rewrite 
them in the p. b.c. form 



sinh(Aj - + ry/2) sinh(A J - + + 27/2) 



N 



sinh(2Aj + 27) 



sinh(A J - - - ry/2) sinh(Aj + - ry/2)J sinh(2Aj - ry) 

sinh(Aj - Afc + ry) 



- n 



-M+l 



sinh(Aj — Afc — 27) 



(A.l) 



where the numbers A-m+i, A-M+2, ■ ■ ■ , Am are related to the v— roots by 



A„ 



-A w , j = l,2,...,M 



(A.2) 



and Revj > 0, thanks to the symmetries of the BAE (2.4). Let us concentrate our 
attention on those BAE solutions which are mostly real, that is those which contain 
an arbitrary but fixed number of complex pairs interspersed in a sea of order iV 
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of real roots. The reason for this restriction will become clear in the sequel. The 
counting function associated to a given solution v\, . . . , vm is defined to be 

Z N (x;vi,..., v M ) =<P-y/2( x + @) + 07/2^ - B) + N'^px) 

1 ^ (A.3) 

j=-M+l 

where 

, . . sinh(io; + x) , . 

4> a {x) = i log — £ (A.4) 

sinn(«a — 

The logarithmic branch in eq. (A.4) is chosen such that 4> a (x), and as direct 
consequence Z^[x ; v\, . . . , vm), are °dd. 

The BAE (A.l) can now be written in compact form 

Z N (\j ;v 1 ,...,v M ) = SttA^ 1 Ij, j = -M + 1, -M = 2, . . . , M (A.5) 

where the quantun integers Ij entirely fix the specific BAE solution, and therefore 
the BA eigenstate, and by construction satisfy = —Ij for j = 1,2, ...,M. 
Notice also that by definition Z^{Q ; v\, . . . , vm)— 0, but A = is not a root, due 
to (A. 2). Rather, A = is always a hole, from the p.b.c. point of view. 

To begin, consider the case when is even, M = N/2 (i.e. J = 0) and all 
roots are real. This state (the groud state for even N) is unambiguously identified 
by the quantum integers 

Ij=j .?' = 1,2,..., TV/2 (A.6) 

and the corresponding counting function is indeed monotonically increasing on the 
real axis, justifying its name. To our knowledge, the existence itself of this BAE 
solution has not been proven in a rigorous analytical way. But it is very easy to 
obtain it numerically for values of A^ in the thousands and precisions of order 10~ 15 
on any common workstation. 
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Next consider removing J roots from the ground state. J is the quantum spin 
of the corresponding new BA state. For 7 sufficently small, one now finds that 

N/2 + J < Z N (+oo;v h ...,v M ) <N/2 + J + l (A.7) 

so that, together with the actual roots satisfying eq. (A. 5), there must exist positive 
numbers x\, X2, • • • , x v , with v > 2 J, satisfying 

Z N {x h ;v h ...,v M ) =2irN~ 1 I h , h = 1, 2, . . . , v (A.8) 

where the Ih are positive integers. If Zn(x) is monotonically increasing, then 
necessarily v = 2 J and the integers . . . , In/2-j-, h, ■ ■ ■ , hj} are all distinct. 
The numbers xi,X2, ■ ■ ■ ,x v are naturally called holes. For J held fixed as N be- 
comes larger and larger, it is natural to expect that the counting function is indeed 
monotonically increasing, and numerical calculations confirm this expectation. For 
larger values of 7 the situation becomes more involved. Numerical studies show 
that, first of all, Z n (x) develops a local maximum beyond the largest root, while 
still satisfying the bounds (A.7), as 7 exceeds a certain (J— dependent) value. 
For even larger values of 7, the asymptotic value of Zn(x) becomes smaller than 
Z* = 2n(l/2 + J/N), but its maximum stays larger than Z* , provided there is 
indeed a root corresponding to N /2 + J (i.e. 1m = N/2 + J). Up to now, vm is 
obviously located where Z^(x) reaches N/2 + J from below, and one could say 
that there exists an extra hole X2J+1 further beyond, where Zjq(x) reaches Z* from 
above. When 7 reaches a certain critical value, the local maximum lowers till Z* . 
At this point the root and the extra hole exchange their places, and for sligtly larger 
values of 7 the hole is located where Zn(x) reaches Z* from below, while the root 
lays further away, where Zn(x) reaches Z* from above. As the numerical calcula- 
tions show, however, this extra hole with the same quantum integer Im = N/2 + J 
of the largest root hole is spurious, since no energy increase is really associated 
to its presence. When none of N/2 — J root has N/2 + J as quantum integer, 
then Zn(x) does never reach Z* for sufficently large 7, and, strictly speaking there 
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are only 2 J — 1 holes. This time, however, we find that the energy increases with 
respect to the ground state in the same way as if there was "a hole at infinity", 
that is a hole beyond the largest root. Thus v can always be regarded to be even, 
when N is even. 

As an important example consider the definite choice J = 1. Then for 7 < 7r/6 
we find Z N (+oo) > Z*, and there are two holes, with 1 < h < h < N/2 + 1. 
Assume h < N/2 + 1 and consider the interval tc/Q < 7 < n/A. The counting 
function has a maximum Z max (larger than Z*) situated to the right of the largest 
root fjv/i-i, as long as 7 is smaller than the critical value 7* at which the maximum 
lowers to Z*. For 7 > 7*, the maximum is still larger than Z* but is located 
to the left of i>jv/i-i- I* 1 an y case, for tt/6 < 7 < 7r/4 the asymptotic value 
Zn(+oo) is smaller than Z* . When 7 — > (n/4)~, then the largest root as well 
as the maximum Z max are pushed to infinity, and Zn(x) is once again monotonic 
with Zn(+oo) = Z*. As 7 exceeds 7r/4, the last root i>_/v/i-i passes, through the 
point at infinity, from the real line to the line with Imv = tt/2. This pictures 
generalizes to arbitrary J with the two special values 7 = 71/6 and 7 = 7r/4 
replaced, respectively, by 7 = it/ (4 J + 2) and 7 = re/ (2 J + 2). In fig. 7 the salient 
portion of the numerically calculated counting function is plotted for J = 1, A = 64 
and a specific choice of I\, 12- In this case we approximatively find 7* ~ 0.2i7r. 

Finally, consider a BAE solution containing, in addition to a number of order 
N of real roots, also a certain configuration of complex roots. In the v— space, these 
complex roots appear either in complex conjugate pairs or with fixed imaginary 
part equal to irr/2, so that the counting function is real analytic: Zn(x) = Z^{x). 
Moreover, it is fairly easy to show, by looking at the value of the counting function 
a real infinity, that the presence of complex roots implies the existence of holes in 
the sea of real roots. For our next porposes, we shall now consider 7/7? irrational, 
so that all v— roots are finite. Denoting with u q , q = 1, 2, . . . , M c the values of the 
complex roots and with M r = M — M c the number of real roots, we now write the 
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derivative of the counting function as 



27rp N (x) = Z' N (x) = F* e (x) + N- 1 [F* (x) + F' c {x) - F' h {x)} + (K * p 6 )(x) (A.9) 



where 



F e (x) =0 7 / 2 Or + e) + 7/2 (x - e) 

F (x) =^(2x) 



v 

F h (x) =0 7 (x) + ^ \-^ X ~ Xn ) + 07 ( x + X n)\ 
n=l 

M r v 

Nps(x) = ^2 $( x - + ^i( x ) + ^2\5(x - x h ) + 5(x + x h )\ 

j=l n=l 



and K* is the convolution defined by 



+oo 



(A.10) 



(K*f)(x)= J dy4>' 1 {x-y)f{y) (A.ll) 

oo 

The so-called "density approach" consists in replacing, as iV — > oo, both pjy an d 
ps with the same smooth function p, representing the density of roots and holes on 
the real line. This function is therefore the unique solution of the linear integral 
equation (eft. eq. (A.9)) 

2np = FQ + N~ 1 (F + F c -F h )+K*p (A.12) 

which can be easily solved by Fourier transformation. Combining this equation 
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with eq. (A. 9), we now obtain, after some simple manipulations 



Pn = p + G * (pn - ps) (A.13) 
where G* = (2tt + K)^ 1 * K* stands for the convolution with kernel 

G ( x ) = f dk c ikx sinh(7r/2- 7 )A: 
K J J 2tt sinh(7r - 7 )A;/2 cosh 7 A;/2 v ' ; 

Finally, a simple application of the residue theorem to the analytic function pjv(l — 
e -iNZ N ^-i pj ug an i n t e g ra tion by parts lead to the following formal nonlinear 
integral equation for the counting function 

Z N = Z + G* L N (A.15) 

where 

1 _ e -iNZ N (x-iO) 



I _ e iNZ N (x+iO) 

L N (x) = -tN- 1 log - -^-t—^ (A.16) 



and Z is the odd primitive of 2np, namely 

Z = Z 00 + N- 1 (2n + K) * (F + F c - F h ) (A.17) 



with 



x 



((2tt + K) * Fq) (x) = 2arctan (^ smh W^ ( A . 18 ) 

\COSh7TC)/7 / 



Eq. (A.15) is a formal integral equation since the knowledge of the exact position 
of holes and complex roots is required in Z. It becomes a true integral equation for 
the ground state counting function. In any case it is an exact expression satisfied 
by Zn where the number of site enters only in an explicit, parametric way, 
except for the positions of the holes and of the complex roots. After having fixed 
the corresponding quantum integers, these parameters retain an implicit, mild 
dependence on N, for large N. 
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We shall now show that eq. (A. 15) is very effective for establishing the result 
(3.6), which is of crucial importance for the calculation of the S— matrix. It is 
sufficcnt to check that the second nonlinear term in the r.h.s. of eq. (A. 15) is 
indeed of higher order in iV -1 relative to the first. To this purpose observe that 
the integration contour of the convolution in eq. (A.15)can be deformed away 
from the upper and lower edges of the real axis, since for sufficently large N no 
complex roots can appear in the whole strip |Imx| < 7/2 and G(x) is analytic 
there. Indeed Z^ tends to as N — > 00, and one can eplicitly check that the 
imaginary part of is positive definite for < Imi < 7/2 and negative definite 
for > Imi > —7/2. This also implies that the contribution to the convolution 
integral is exponentially small in N for all values of the integration variable (let's 
call it y) where ImZoo of order 1. For |Rey| of order log TV, we find ImZoo(y) of 
order iV -1 , so that the nonlinearity Ljv, rather than exponentially small, is also of 
order N^ 1 . But now the exponential damping in y of the kernel G(x—y) guarantees 
that the convolution integral is globally of order N~ 2 or smaller, provided \x\ is 
kept smaller than (7/vr) log N. Hence we can write 



Finally, the coefficent of the iV -1 term of Z, in eq. (A. 17), can be calculated in 
the N — > 00 limit, with the techniques described at length in paper I. After some 
straightforward albeit cumbersome algebra, this yields 



lim ((2n + K)*(F + F c -F h ))(x) = F(x;x 1 ,...,x M ;Xi,---,qX M ) ( A - 2 0) 



where the higher-level quantity F is defined in eq. (3.8). Together with eqs. (A. 19) 
and (A. 17), this proves eq. (3.6) of section 3, as claimed. 

Let us now consider the problem of calculating the energy of a given BA state, 
through eq. (2.10). We rewrite first the "bare" energy function eo(x) as 



Z N = Z + 0(N- 2 ) 



(A.19) 



eo(x) = -2-k + (f> 7 / 2 (x + 0) - 7 / 2 Or - 0) 



(A.21) 
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Then we calculate 



M +?° M c v 

^0 7/2 (a;) =\N j p s (x)<i> l/2 (x) + 5^0 7 / 2 K) - Xl^/zW - ^7/2(°) 

J = l 4q «=1 fc=l 

+?° M c 

-00 5=1 /i=l 

+00 

+±iV | - pjv(x)] ((1 - G) * 7/2 ) (s) 

—00 

(A.22) 

Through the residue theorem and an integration by parts, the last term can be 
transformed, as done before for the counting function, into an integral of the non- 
linear term Ljy, namely the integral 



+00 

_ AT 

2 



i N I l n( x > 



7 cosh nx fry 



By the same argument used above, this last expression is globally of order TV -1 . 
Finally, inserting the explicit form of the continuum density p(x) into eq. (A.22) 
and recalling eqs. (2.10) and (A. 21), for the energy we obtain eq. (3.1) of the main 
text 

V 

E = Eq + a' 1 X e ( x h) + 0{a- l N- 1 ) (A.23) 
h=i 

where 



E = -na 1 N+a~ 



00 

, ^ - sinh(7r — 37)fc/4 
2N cos kO + 1 + 



1 f dk smh(Tr — 7)fc/2 sin/c6 







k sinh7r/c/2 cosli7/c/2 



sinh(7r — 7)fc/4 
(A.24) 

is the energy of the ground state. 

We would like to close this appendix with a comment on the limitations of the 
density approach, where one deals only with the solution p of the linear equation 
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(A. 12). Regarding p(x) as the actual density of real roots and holes in the N — > oo 
limit, it is natural to use it to replace summations with integrals. What one learns 
from the exact treatment presented above as well as from computer calculations, 
is that the error made in such a replacement depends crucially on the large x 
behaviour of the quantity which is to be summed. This error is down by N^ 1 only 
when there is exponential damping in x. This means, for instance, that the integral 
of p(x) does not reproduce in general the exact number of real roots and holes, but 
rather some 7— dependent quantity close to it. Misundertanding this for the actual 
number of real roots and holes would lead to the absurd result that the holes have 
a 7— dependent value of the SU(2) q spin, which is instead necessarily integer or 
half-integer and, in the particular case of the holes, just 1/2 for any value of 7. 
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8. Figure Captions 



Fig.l. The "ground state dominating" configuration of local height variables. 

Fig.2. Graphical representation of the vertex-face correspondence. The strip sand- 
wiched between the two time lines represents the unit time evolution. 

Fig.?,. Graphical representation of the one-soliton operator corresponding to the 
particle with rapidity 82 in a system with v — 4 particles. At each intersec- 
tion, the appropriate two-body S— matrix acts. In the collisions against the 
wall, 62 is flipped and the two-by-two matrix ^(i^) acts on the internal 
states of the soliton. The choice of 62 as largest rapidity is done purely for 
graphical convenience. 

Fig A. One of the local height configurations that dominate the one-hole BA state, 
of spin J = 1. 

Fig.h. The "renormalized" version of Fig.4. 

Fig.6. Flipping the last portion of the configuration from the solid to the dotted 
line does not change the quantum spin J = p/2 — 1, but tranform the type 
II state into type I. 

Fig.l. Plot of the N/2tt times the counting function Zjy(x) versus tanh x for N = 64, 
J = 1, 6 = .15 and various values of 7 from n/Q to 0.9997r/4. The quantum 
integers associates to the two holes are I\ — 8 and I\ = 21. The critical value 
of (N/2n)Z N is N/2 + J = 33, while that of 7 is, roughly, 7* = 0.21tt. 
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